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Abstract. We study the splitting probabilities for a one-dimensional Brownian 
motion in a cage whose two boundaries move at constant speeds ci and C 2 . This 
configuration corresponds to the capture of a diffusing, but skittish lamb, with an 
approaching shepherd on the left and a precipice on the right. We derive compact 
expressions for these splitting probabilities when the cage is expanding. We also 
obtain the time-dependent first-passage probability to the left boundary, as well as the 
splitting probability to this boundary, when the cage is either expanding or contracting. 
The boundary motions have a non-trivial impact on the splitting probabilities, leading 
to multiple regimes of behavior that depend on the expansion or contraction speed of 
the cage. In particular, the probability to capture the lamb is maximized when the 
shepherd moves at a non-zero optimal speed if the initial lamb position and the ratio 
between the two boundary speeds satisfy certain conditions. 
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1. Introduction 

A lamb escapes from a farm and has the bad idea to roam near a precipice. The shepherd 
wonders what is the best strategy to catch her lamb alive. Indeed, the lamb wanders 
randomly when the shepherd stays still. However the lamb is skittish and moves away 
from the shepherd and toward the precipice whenever the shepherd approaches. Shonld 
the shepherd stay still and hope the lamb will come to her, or shonld she walk toward 
the lamb and hope that she reaches the lamb before it goes over the precipice? 

ac c Cj C2 

shepherd lamb precipice 

(a) (b) 

Figure 1. Equivalent formulations of the lamb capture process, (a) Shepherd’s 
reference frame: the lamb and the precipice drift toward the shepherd at speeds ac and 
c. (b) Reference frame where the lamb only diffuses: the shepherd and the precipice 
approach the lamb at speeds ci = ac and C 2 = (1 — a)c. 

If the shepherd walks toward the lamb with speed c, we assnme that the lamb 
moves away at a slower speed (1 — a)c, with 0 ^ a ^ 1, that is snperimposed on its 
diffnsive motion. In the reference frame of the shepherd (Fig. [^a)), the lamb diffnses 
and approaches the shepherd with speed ac. However, the precipice also approaches the 
lamb with speed (1 — a)c. Initially, the lamb is at Xq, the precipice is at Lq, while the 
shepherd is hxed at the origin. The probability that the lamb and shepherd meet before 
the lamb goes over the precipice coincides with the splitting probability for the lamb to 
eventnally reach the origin. In tnrn, this latter problem is eqnivalent to the splitting 
probability for nnbiased diffnsion that starts at x'q = xq —To/2 to reach the left edge of 
an asymmetrically contracting cage, whose left and right edges are at —Lo/2 + cit and 
To/2—C 2 t, respectively, with ci = ac and C 2 = (1—a)c (Fig. [^b)). While the contracting 
cage is the relevant sitnation for the lamb captnre problem, we also investigate the case 
of the expanding cage. Althongh the valne of a does not need to be restricted, we stndy 
the range a G [0,1], where the splitting probability has the richest behavior. 

Conditional exit from a hxed interval is a classic problem of random-walk theory 
The role of a moving bonndary has been considered recently, e.g., a diffnsing or 
an oscillating trap at the edge of the interval mil- For ballistically moving bonndaries, 
which is onr focns, the inhnite-time snrvival probability in an asymmetrically expanding 
cage and the time-dependent snrvival probability in a symmetrically expanding cage 
were investigated by Bray and Smith j^|^. Here we derive exact expressions for the 
hrst-passage probability at any time and the splitting probabilities to one of the edges 
of the cage when each wall moves ballistically at an arbitrary speed. Contracting and 
expanding cages are considered. The snrprising conseqnence of onr resnlts is that the 
splitting probabilities depend non-monotonically on the speed c for a range of initial 
conditions and a valnes. 
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In Sec. we derive the splitting probability in an expanding cage. In Sec. we 
determine the hrst-passage probability at any time and the splitting probability to one 
of the boundaries for both an expanding and contracting cage. Finally, in Sec. we 
demonstrate that the splitting probability can have a non-monotonic dependence on 
the speed c. Using these results, we answer the shepherd’s question of what is the best 
strategy to capture the lamb. 


2. Expanding Cage 


It is convenient to treat the problem in the reference frame of the left boundary 
(Fig. 13 a)), where the Brownian particle drifts to the right with speed ac and the 
right boundary also drifts to the right with speed c. Thus the right boundary is located 
at L{t) = Lo+ct > Lq. We focus on the splitting probability to reach the left edge of an 
expanding cage Lq) as a function of the initial particle position xq and the initial 

interval length Lq (with TZ^ the splitting probability to the right edge). Following the 
approach of Ref. j^, the backward Fokker-Planck equation for the splitting probability 
is 


D 


dxl 


dC^ dC^ 


ac- 


= 0 , 


dxo OLq 

or, by introducing the rescaled variables y = cxq/D and A = cLq/D, 


(la) 


da da 

-h O. -1- 

dy^ dy d\ 


(lb) 


with 0 ^ 1 / ^ A and the boundary conditions £®(0, A) = 1 and £®(A, A) = 0. 
In the spirit of j^, we seek a solution of the form 


a{y, A) = ^ 


( 2 ) 


The left boundary condition U®(0, A) = 1 leads to 

U^(0, A) = + 6n)e-"("+“)" = 1, (3a) 


while the right boundary condition £®(A, A) = 0 gives 

a{\ A) = ^ = 0 . 

nEZ 

By shifting the index of the second sum, n —>■ n — 1, and after some simple algebra, the 
condition above can be written as 


r (A, A) = = 0 . 

nEZ 


(3b) 
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Since Eqs. ^ hold for all A, we obtain the following relations for the coefficients in ([^; 
flO + ^0 = 1 ! 

an + bn = 0 Vn 7 ^ 0, (4) 

On + bn-1 = 0 Vn. 

If the initial cage length Lq —)■ oo and the initial position of the particle is far from 
either boundary, then the splitting probability —>■ 0. Using these facts in ([^ imposes 
Oq = 0. Together with the relations Q, we hnally obtain the splitting probability to 
the left edge: 


C{y, A) = + ^ [« 


-{n+a)y _ nyl -n{n+a)\ 


( 5 ) 


n=l 


Note that the splitting probability to the right edge can be obtained from g by 

K'(!/.A)1„ = £=(A-!/.A)|,_„. 

To obtain the solution in the second formulation where an unbiased Brownian 
particle starts at Xq = xq — Lq/2 in a cage whose boundaries are located at —Lo/2 — cU 
and Lo/2 + C 2 f at time t (Fig. [^b)), we replace a, c and Xq with their corresponding 
expressions in terms of Ci, C 2 and Xq to give 


i:"(x'o,Lo) = e-"i(2x'+Lo)/2D ^ E g-n[n(ci+C2)+ci]Lo/D 

n=l 

X < p[-^(ci+C2)+ci](2a;^+Lo)/2X' 


[n(ci+C2)](2a;Q+Lo)/2D 


( 6 ) 


When the initial interval length Lq 3> 1, the splitting probability is well approximated 
by its hrst term 

To) ~ ^ ^7^ 

which is exponentially larger than all higher-order terms in the series in ([^. 

As a useful counterpoint, we can recover this last result by applying the “free 
approximation” [4, 10 , in which the concentration within the interval is assumed to 
retain the same Gaussian form as a diffusing particle on the inhnite line with no imposed 
boundary conditions. This approximation relies on the boundaries being outside the 
range where the probability distribution is appreciable. Thus we assume that the 
concentration prohle is 

v47riAf 

where the unknown amplitude A{t) accounts for the loss of probability within the domain 
and should be determined self consistently. In this free approximation, the flux at any 
point in space is 

. = g-(x-x')V4Di 

dx y/WnDO 

Thus the total flux leaving the cage is the sum of the flux at the two boundaries: 
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dA 

-r = -^ 
dt 




Integrating this equation to finite time, the amplitude is 

^—x'r, — C2t\ g-ci(-Lo+2a;^)/2n> 


( 8 ) 


p—C2{Lq—2x'q)I2D 

In A{t) = -erfc I 


^/Wi 


erfc 


^+a;o-Cif 

^/ADi 


(9) 


For f —)■ cx), this expression reduces to 


A(t = oo) = exp 


^-C2(Lo-2x'q)I 2D _ -cPLo+2x'g)/2D 


1 _ ^-C2{Lo-2x'q)/ 2D _ ^-cPLo+2x'g)/2D 


Ln —y oo. 


( 10 ) 


Here A{t = oo) represents the large-time limit of the survival probability, namely, the 
probability that the particle has not touched either of the boundaries. It can thus be 
written as H(f = cx)) = 1 — C^{xq, Lq) — 71'^{xq, Lq), which Anally leads to 


C^x'q, Lo) ~ er^dU+‘lx',)l2D _ ^-C2{Lo-2x',)I2D 


( 11 ) 


in agreement with Eq. Q. 

Unfortunately, this backward Fokker-Planck approach does not seem to be 
adaptable to a contracting cage. In this case, the general solution involves a complex 
exponential dependence in y. Therefore, the device used to simplify the right boundary 
condition to the form given in (|^ no longer holds. To obtain the splitting probability 
in this case, we apply a more general framework following the methods of [9,10 . 


3. Contracting and Expanding Cage 


We now turn to the general case of a cage that can be either contracting or expanding. 
Again, we study the problem in the reference frame of the left boundary that is fixed 
at a; = 0. For the expanding cage, the right boundary moves to the right at speed c 
and the particle drifts to the right with speed ac in addition to its diffusion. For the 
contracting cage, the right boundary and the particle both drift to the left. Let x denote 
the position of the particle at time t. We first compute the propagator p(x, t) for the 
particle in this cage, viz., the probability for the particle to be at x at time f, by solving 
the forward Fokker-Planck equation (see, e.g., |^) 


dp dp d‘^p 

dt dx dx'^ ’ 


( 12 ) 


with the initial and boundary conditions p{x, 0) = 5{x—Xq) and p(0, t) = p(^L(t),t) = 0, 
and with L{t) = Lq ± ct. Here the upper sign corresponds to the expanding cage and 
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the lower sign to the contracting cage thronghout this section. For the contracting cage 
the process stops when the two bonndaries meet at t = Lq/c. 

When the boundaries are immobile, so that the cage length L{t) = Lq is constant, 
the elemental solutions of the Fokker-Planck equation with these boundary conditions 
are well known 


fn{,x,t) = sin 


nnx 


exp 


acx a^c^t n^n'^Dt 


2D 


AD 


LI 


with n eN. 


To account for the interval length changing linearly with time, we follow the method 


employed in and adapted by Bray and Smith [9] and postulate a solution to Eq. (12) 
with absorbing boundary conditions at x = 0 and x = L{t) of the form 

acx a‘^cH\ ( 2 


Pn{xO) = fi'(x,f)sin 


( nvrx \ 


exp ± 


2D 


AD 


\L{t)J 

Substituting this trial function into Eq. (12), we obtain 

/ mix 


exp - n 71 D 


Lm 


d'^g dg 


nil 


L{t) V dx L{tf 


(13) 


(14) 


We notice, as done in 0, that we can seek a form for g{x,t) that cancels both the 


left- and right-hand sides of (14). Thus g{x,t) niust simultaneously solve 

d^g ^ 

dt ’ 


dx‘^ 


^„dg cx 
2D^ = =F—-0 . 
dx L{ty 

These equations imply that the function g{x,t) has the form 

K 


g{xO) = 


Vm 




(15) 


with K a constant. The general solution can now be written as a superposition of the 
basis functions p„(x, t)\ 

a^c^t n^TT^Dt' 


P{x,t) = 


n£N 


Ctn 

7W) 


sm 


mix 

Lit) 


exp =F 


9 

X C 


ADLit) 2D 


acx 

± -rw I exp 


AD 


LoLit) 


(16) 


For the initial condition p(x, 0) = S{x — xq), we use the identity 


E 

neN 


. f n7ix\ . f n7iXo\ Lq 


for 0 ^ X, Xq ^ To to hnally obtain the solution for the given initial condition as 

2 . / mix\ , / n7rxo\ 


= Y 


sm 


nSN 


-y/LoT(f) \L{t) J \ Lq J 

c(x^ — Xq) ac(x—Xo) a'^cH n^7i‘^Dt\ 

ADLit) 2/1 AD~ ~ LoLit) ) 


X exp =F 


(17) 
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F(0,i) = C^[ , 

OX la;=0 




2n7TD 

^/LJm 


sm 


n7iXo\ 


f 


± 


CXn 


V 4DLc 


T 


acxo 

JlD 


a^cH 

AD 


LoL{t) J 


(18) 


The splitting probability to the left edge is the time integral of this first-passage 
probability. As noted previously, the temporal integration range depends on the sign of 
the speed. Finally, the splitting probabilities for the contracting and expanding cage, 
and respectively, are 


£'^(xo, Lq) 



F(0,t) 


E 

neN 


2n7rD 
—-;= sin 
cy/Lo 


rnixo 

Lo 


^-c{xo-aLo)'^/ADL q 


^ ^DDD/cLo 



exp 



n^7r^D\ 
cL )’ 


(19a) 


r(x„.L„) = /■“* F(O.t) = 5 : ^sm 


nEN 


^ ^-DDd/cLo 


V Lq / 
dL 


' Lq 


L3/2 


exp ( — 


a'^cL 

+ 


AD 


cL 


(19b) 


Expression (19b) for the splitting probability in an expanding cage can be shown 
numerically to perfectly match the simpler form ([^ derived by the backward Fokker- 
Planck equation. 

However, these splitting probabilities are not convenient for numerical evaluation. 


Instead, it is expedient to use the Poisson summation formula 11 

^h{n) = ^ h(27rm), 

riEZ mEZ 

with h{x) = dt e“*®^/(f), to give the alternative expression 

I'Lo 


C^{Xq, Lq) — 


mEZ 


^ ^—c{xq—olLo)‘^/ADLo 

AttD 


dL 




^orcLjAD 


X exp 


cL(4Lgm^-l-a;o) 

ADLq{Lq-L) 


xq cosh 


cLxom 

D{Lo-L) 


— 2mLn sinh 


cLxom 

D{Lo-L) 


. ( 20 ) 


which is more suitable for numerical evaluation (and similarly for C^{xo, Lq)). 


4. Optimal Capture Criterion 

We now turn to our original question; what is the optimal strategy for the shepherd to 
catch her skittish lamb without driving it over the precipice? In the shepherd’s reference 
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frame, the lamb approaches at speed ac while the precipice approaches at a higher speed 
c. The probability to catch the lamb in this contracting cage is the splitting probability 


C^{xo,Lo) (19a). What speed c maximizes this splitting probability? 

Partial conclusions can be drawn by studying the limits of c —)■ 0 and c —)■ oo. If 
c = 0, the splitting probability £'^(xo,To) is a linear function of its initial position (see, 
e-g., 0) 


C^{xo, Lq) — 


Lq — Xq 


The qualitative behavior of C^{xq,Lq) when c —)■ oo can also be easily understood. In 
this limit, if the time xq/{ ac) for the lamb to reach the shepherd is smaller than the time 
{Lq — xq)/[{I —a)c\ for the precipice to catch up to the lamb, then £^(xo,To) 1, while 

C^{xo, Lq) —)■ 0 otherwise. These two times match when a = xq/L. Thus in the limit of 
large speed, the splitting probability reduces to a step function, with O'{xo, Lq) ^ 1 for 
Xq/L < a and £^(xo, Lq) 0 for Xq/L > a . 

Let us now focus on the hrst order of the splitting probability for small speeds. In 


this case, the integral over L in Eq. (19a) 


rLo 


1 = 


dL 

L^ 


exp 


a^cL 


AD 


cL 


can be developed with respect to c by expanding exp{a^cL/AD) and integrating by 
parts. This yields 


exp 


Then using 


cLq 


E 

n=l 

oo 


,2r3/2 


c^L 






0 




( 21 ) 


-1)^ 


sm nz = —- 


n 


z 

2 ’ 


-1)^ 


n=l 


sm nz = 




z 

12 


Tl'^Z 

u 


we obtain the small-speed behavior of the splitting probability: 

Lo 


C^’%Xq,Lq) = 


Xq cxq {Lq - Xo) [(So; - l)Lo - . . 

T ttvt- 79 -^ o { c ). 


LI 


( 22 ) 


Lq 6D 

Here we also quote the limiting form splitting probability for the expanding cage, which 


is found from (19b) by the same steps as outlined above. 


As a result of the a and xq dependence of the first-order term in c given above, 
C^{xq,Lo) is an increasing function of speed at c = 0 when a > ^(^Xq/Lq + 1) and 
decreasing otherwise. Combining this fact with the limiting behavior for c —>■ oo, we 
deduce that the splitting probability can be a non-monotonic function of the speed, for 



















Optimal strategy to capture a skittish lamb wandering near a precipice 


9 



Figure 2. Splitting probability phase diagram for the (a) contracting and (b) 
expanding cage in the speed ratio (a) and initial position (xq) plane. In (a), « 1 

for c —>■ oo above the solid line, while « 0 for c —>■ c» below. Above the dashed line 
is an increasing function of c at c = 0 and a decreasing function below. These lines 
delineate four zones of behavior as discussed in the text. In (b), there are two zones. 


specific values of Xq and a. This leads to rich behaviors for the splitting probability, as 
illustrated in Fig. |^a). 

We can now give advice to the shepherd. There are four distinct strategies, 
corresponding to the four zones of Fig. |^a): 


(i) “Dangerous zone” (lower right). Here, either the lamb is very fearful, a 1, or 
sufficiently close to the precipice, a < |(xo/To + l) and a < Xq/Lq, so that 
monotonically decreases with c. Thus the optimum strategy for the shepherd is to 
not move. 


(ii) “Safe zone” (upper left). Here, either the lamb is not very fearful or is sufficiently 
close to the shepherd, a > |(xo/To + l) and a > Xq/Lq, so that monotonically 
increases with c. The shepherd should run as fast as possible to maximize the 
probability to catch the lamb. 

(hi) “Optimizable zone” (upper right). Here, the lamb is close to the precipice but not 
too fearful, a > |(xo/To + l) and a < Xq/Lq so that has a maximum as a 
function of c. Thus there is an optimal speed that maximizes the probability for 
the shepherd to catch the lamb. 

(iv) “Dilemma zone” (lower left). Here the lamb is close to the shepherd but also very 
fearful, a < ^{xq/Lq + 1) and a > Xq/Lq. Thus initially decreases with c before 
eventually increasing. Thus if the shepherd is unht, she should stay still. However, 
if she is sufficiently £t, she should run as fast as possible. 


For the expanding cage, the splitting probability given in Eq. (22) is now an 
increasing function of speed at c = 0 for a < \{xq/Lq + 1). Using this small-speed 
dependence, together with the limiting behavior C^{xq, Lq) —)■ 0 for large speed, gives the 
phase diagram shown Fig. [^b) . There again exists a zone in the phase diagram where the 
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splitting probability can be maximized with respect to the speed. Here, the maximum 
of the splitting probability arises from the interplay between two competing effects from 
the cage expansion. Indeed, consider the complementary probability that the lamb 
never reaches the left boundary. The larger the speed, the lower the probability for the 
lamb to reach the left boundary. At the same time, the larger the speed, the higher the 
probability for the lamb to also not reach the right boundary, which implicitly increases 
the probability to reach the left boundary. As a result of these competing effects, there 
exists an optimal speed of expansion that maximizes the splitting probability. 




Speed c 

(b) 


Figure 3. Splitting probability to catch the lamb in the (a) “optimizable zone” and 
(b) “dilemma zone”. In (a), the splitting probability is increased substantially (by more 
than a factor 3 for a = 0.92, Xq = 19, and Lq = 20) when the shepherd moves at the 
optimal speed Copt instead of staying still. In (b), the probability to catch the lamb is 
diminished if the shepherd runs slower than speed c*. Here the parameter values are 
a = 0.08, a;o = 1 and Lq = 20). 


5. Conclusion 

We analytically determined the splitting probabilities for a one-dimensional Brownian 
motion in a cage whose boundaries move at constant speeds Ci and C 2 . We analyzed 
both the cases of contracting and expanding cages. In addition, we calculated the time- 
dependent hrst-passage probabilities at each of the boundaries. Intriguing behaviors 
of the splitting probabilities arise as a consequence of the ballistic boundary motion. 
Indeed, we found that the splitting probabilities can vary non-monotonically with the 
relative speeds of the boundaries, depending on the initial position of the Brownian 
particle and the ratio between the two boundary speeds. In the context of a fearful lamb 
near a precipice (Fig. [^, this non-monotonicity dehnes a non-trivial optimal strategy 
for the shepherd to catch the lamb. 

This approach could be extended to determine the splitting probability to a subset 
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of a growing d-dimensional sphere for a particle that starts somewhere within the 
interior. It would also be interesting to extend the conditional exit problem to the 
case of non-linear displacements of the boundaries where the lamb capture probability 
should also have a non-trivial optimization. 
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